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"Microwave & Infrared Observations of Molecular 
Spectra in Comet Kohoutek" 

Funds were granted to assist with travel and other expenses 
connected with observations on Comet Kohoutek in both the microwave and 
the infrared regions. 

Microwave measurements were made with a maser amplifier on 
the 80 ft. radio telescope of the Naval Research Laboratory in conjunction 
with personnel there. Nine transitions of possible molecules, ^O, NH^, 
CHgOHj, N^O, and OH were searched for in Comet Kohoutek in the 
frequency range 22.2 - 25.2 GHz. These molecules were not detected, 
but upper limits for the optical depth, mean column density and production 
rate were derived for each of the molecules. In view of the reported 
detection of CH^CH, the upper limits obtained for ammonia, for example, 
were surprisingly low. A detailed report of this work will be found in 
the accompanying manuscript, which has been submitted for publication in 
the journal ICARUS. 

A search was made for the resonant scattering of solar infrared 
light at 5^i by CO molecules in Comet Kohoutek. Two observing periods 
were scheduled on the Lick 120 in. telescope in October and November, 
but unfortunately cloudy weather during both periods prevented any obser- 
vations, Additional observation periods were scheduled on the 84 in. 
telescope at Mauna Kea. Again, there were some weather difficulties 
which, with logistic problems with equipment and cryogenic materials 
limited the observations. However, observations of moderate quality were 
obtained. These set an upper limit on the evolution of GO from Comet 
Kohoutek at about the level of expected rates. It is still quite low, however, 
in comparison with the reported rate of evolution for CH^CN. A detailed . 
description of these results is included in the accompanying manuscript 
which is also scheduled for publication in ICARUS. 
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ABSTRACT 


' ' t 

Nine transitions of the possible parent molecules H^O, 

♦ 

NH q , CH^OH and N 9 0 as well as the OH radical were searched for 
in Comet KohouteK (1973f) in the frequency range 22.2-25.2 GHz. 

V 

These molecules were not detected but the upper limits for the 
optical depth, mean column density and the production rate are 
derived for each of the molecules. These results are discussed 
and compared with the reported detections of HCN and CHgCN 
emission and OH absorption. • ‘V 



The early discovery of the bright, young Comet Kohoutek 
(1973f) provided an excellent opportunity to search for several 
of the possible parent molecules (Huebner, 1970a, 1970b) 
evaporated from the icy nucleus of the comet near perihelion 
passage on December 28. During the period November 1973 - 
January 1974 we conducted a search for nine microwave transitions 
of five molecules — ^O, NH^, OH, CH^OH, and N 2 O* The searches 
yielded negative results. 

The 85-foot (26m) reflector of the U. S. Naval Research 
Laboratory at Maryland Point Observatory was used for all 


observations. For the ammonia, methanol and nitrous oxide 
searches the receiver was equipped with a maser preamplifier 


(T- A ~ 
v System 

Berkeley . 


200°K) developed at the University of California at 
For the water vapor, hydroxyl radical and some of 


the ammonia searches the receiver used a tunable parametric 
amplifier (T„ , 850°K) made available by the National Radio 

Astronomy Observatory.'*' Alternate back ends consisted of a 
bank of fifty 50 kHz contiguous filters and parallel banks of 


twenty-five 100 kHz and twenty 500 kHz filters with the same 
center frequency. At 24 GHz this corresponds to velocity * 
resolutions of approximately 0.6, 1.2 and 6 km/sec respectively. 
Table I contains a summary of the observations. 


^The National Radio Astronomy Observatory is operated by 
Associate Universities, Incorporated, under contract with the 
National Science Foundation. 
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The nine transitions cover the frequency range 22,2-25.2 GHz. 
At these frequencies the antenna has a half-power beamwidth of 
— 2. 2-2,0 arcmin, an aperture efficiency of approximately 0.4 
and a main beam efficiency of approximately 0.5, An on-off 
technique with beam switching between comet and adjacent sky 
was used for all the data with typically 5-10 minute runs on 
source. An ephemeris for the topocentric position and velocity 
of the comet was generated at NRL using periodically updated 
orbital elements provided by B. G. Marsden (Smithsonian- Astro- 
physical Observatory) . The observat ions concentrated on the 
predicted positions for the comet nucleus but some of the 
searches included a grid around the predicted position and the 
extended tail of the comet. 

Assuming the comet to be optically thin with a negligibly 
small optically thick center., the upper limit for the optical 
depth., r, in terms of the observed peak-to-peak antenna 
temperature, T^, is 

t A = 7 ?B n(T " T b ) 0--e~ T > = V B B(T-T b )r (1) 

where pg - main beam efficiency, B = beam dilution factor, - 
2.7°K microwave background continuum and T - temperature of 
cometary gas. A value of T - 200°K was assumed for the evaluation 
of the upper limits. This value is consistent with the 200°K 
equilibrium temperature of a water clathrate model for the comet 
at a heliocentric distance of R = 1 A.U. (Jackson, 1972), 

The optical depths at the center of a Doppler-broadened 
line for the linear molecules and OH, the symmetric top 



NH 3 and the asymmetric tops HgO and CH^OH were taken to be 
respectively (Huebner, 1970b), 


r 
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r 
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16 ir 3 f y B\/C |j 3 v 3 (2J + 1) K 3 _ w / kT ^ 

3ck B/2 T 5 ^ A v J (J + 1) 
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( 2 ) 


where all quantities are in cgs units, h = Planck's constant, 

k ~ Boltzmann's constant, c — light velocity, f y — 1 for 

transitions considered, all in ground vibrational level, 

V = frequency (Hz), jx ~ dipole moment (esu~cm), T » cometary 

gas temperature (°K), A v — line half-width at e ^ of maximum 

intensity (Hz), A,B,C are rotational constants (Hz)-, A = B = C, 

J and K are quantum numbers of upper rotational sate, W = 

energy of upper rotational state, S = line strength of transition 

-2 

and N = column density of molecular species considered in cm 
The beam dilution was taken into account in the usual 
*' way using the optical depth averaged over the antenna beam < r > 

= B r in Eq. (1) and the column density derived from Eq. (2) is 
then the uniformly smeared out or averaged column density of the . 
molecules in all states in the antenna beam ( N ) (Huebner and 
Snyder, 1970) « Table II contains the results for the upper 
limits of the average optical depths and column densities. 

Since the velocity of the possible parent molecules such 
as H 2 0, NH 3 , CH 3 0H and N (> 0 is determined by the equilibrium . . 
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temperature with the icy neculeus (Jackson, 1972) the linewidths 
were calculated from the expression for the Doppler half-width 
at e7^ of maximum intensity A v~v \/"2kN T/M /C where T = 200°K, 

M = molecular weight and = Avogadro's number.' The full-width 
at half-intensity 2 (ln2) 1//2 A v for these molecules would be 
38-58 kHz. The Doppler width was also used for the OH radical 
upper limits, however, the velocity of OH depends upon the 
actual radical production mechanism (Jackson, 1972). OH line 
widths of about 2 km/sec may have to be considered (Turner, 1974). 

Upper limits for the molecule production rate Q in molecules/ 
sec/steradian were estimated from its relation to the average 
column density < N > (Huebner and Snyder, 1970) . 

Q = v - fg - ° 2 [s cos" 1 § - (r® - s s ) w + r )- 1 < N > (3) 

o . 

where v is the average expansion velocity of the molecular species 

which is approximately 0.3-0, 4 kra/sec, - destruction radius 

for the molecule where r = v t with r the lifetime of the 

o o o 

molecules, A = geocentric distance of comet, 8 = full width of 

antenna beam at half power and s is the smaller of the two • 

quantities r Q and A (8/2). Typically, A (0/2) = 4-5 x 10^ km. 

At a heliocentric distance of 1 A.U. the' parent molecules have 

4 5 

lifetimes against phoiodissociation on the order of 10 -10 sec 

(Jackson, 1974) . Values of the molecular yield Q (molecules/sec 

4 5 

/ster) estimated for the typical values r = 10 and 10 km are 

o* 

contained in Table III. 

An estimate for the molecular yield of OH can be made 
if one considers that the OH radical is in an extended cloud 
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around the nucleus with an r ~ 10 km at H = 1 A.U. (Delsemme, 

* o 

5 

1973), or r 0 ~ 6 x 10 km at R = 0.8 A.U. appropriate for our 

w ° m 

OH observations. In this case the OH filled the beam and the 

16 -2 

upper limit of ( N ) < 4.2 X 10 cm yields a value of 
31 

Q < 2 x 10 OH radicals produced/sec/steradian. Turner 

(1974) has reported the detection of the 1665 and 1667 MHz lines 

of OH in Comet Kohoutek in absorption against a 3°K background, 

T4 _9 

He derived a value of N = 6 x 10 cm w in the ground state with 

a beamwidth of 18 arcmin and estimates the size of the OH halo 

to be about 11 ai'cmin in early December which corresponds to 
5 

r o ~ 4 x 10 km in agreement with the value used above. 

Huebner, Snyder and Buhl (1974) have detected HCN in Comet 

Kohoutek with estimates for the HCN production rate of 
26 

Q ~ 2-7 x 10 molecules/sec/steradian. If the production 
rates for the molecules of Table III are actually similar to 
that of HCN then our upper limits for Q are not very stringent 
* except for the case of. NH^ whose production rate might be 
comparable to HCN. Our non-detection of NH^ may have been 
** due to a smaller, destruction radius or to a relatively lower 
initial abundance within the icy nucleus. The contrast is 
more striking when one compares NII^ with CH^CN. Ulich and 
Conklin (1974) have reported the detection of two vibrationally 
excited rotational transitions of methyl cyanide, CII^CN in 

16 

Comet Kohoutek with an estimated column density of N = 3 x 10 
—2 

cm in early December. For the beam dilution used (B = .06) 
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15 ♦ —2 

this would correspond to ( N ) = BN = 2 x 10 cm . For 

4' 

typical values of the,, destruction radius, = 10 km, this 

* 29 

would be equivalent to a production rate of Q = 2 x 10 

5 

CH 0 CN molecules/sec/s teradian while for r - 10 1cm, the 
j o 

28 

rate would be Q = 4 x 10 . Our upper limit indicates that 

CHgOH was produced at a slower rate than the detected CH 3 CN 
but the production rate of CH^OH could have been comparable to 
that of the detected HCN. 

The temperatures required to produce 1 mm of gas pressure 
for N 2 0, NH 3 , HCN, CH^CN, CH 3 0H and HgO are respectively ' 130° , 
164°, 202° , 226°, 229° and 266°K (Handbook of Chemistry and 
Physics, 1967-1968) * The two molecules detected are among the 
least volatile on this list. Our upper limits set for the 
abundances of N 2 0 ' CH 3 0H and HgO are sufficiently large to 
perhaps not be very significant. Indeed, the presence of 
H 0 0 in the cometary gas is inferred by Wehinger et al. (1974) 
who detected the H 2 0^ ion in the comet tail. However, the low 
abundance of NH 3 is surprising and may be due to its having 
been lost by evaporation during some past history of the comet 
or to formation of the comet at temperatures which were too high 
to condense much NH 3 in the nucleus. 

We gratefully acknowledge B. G. Marsden for his many updates 
of the comet orbital elements. We thank P. K. Schwartz, J. M. 
Bologna and T. P. McCullough for their assistance with the 
receiver' sys terns and observations, "W. M. Jackson for his advice 
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and W. F. puebner, L.E, Snyder and D. Buhl for making available 

their HCN results before publication. We are especially grateful 

# 

to Professor C. II. Townes for his helpful discussions. 
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TABLE I 


MOLECULE SEARCHES IK COMET KOHOTOEK (1973f) 


Koleevls 

Transition 

Rest 

Frequency 

KHz 

Observation 

Dates 

1973-1974 

A(AU) 1 

Filter 

Resolution 

kHz 

Receiver* 2 

Integration 

Time 

(minutes) 

H z O 

6 16~ 5 23 

22235.08 

Nov 14, Jan 2,18 

0.9 

50 

P . 

* " 120 

H 2° 

6 16" 5 23 

22235.08 

Nov 14, 16, 26 

1.6 

100,500 

P 

140 

nh 3 

(1,1) 

23694.49 

Dec 22 

1.1 

50 

M 

85 

nh 3 

(1,1) 

23694.49 

Nov 23 

1.5 

100,500 

‘ P 

90 

nh 3 

(2,2) 

23722.63 

Nov 23 

1.5 

100,500 

P 

20 

nh 3 

(3,3) 

23870.13 

Dec 16, 20, 23 

1.0 

50 

/ 

M 

140 

nh 3 

(3,3) 

23870 , 13 

Nov 19 

1.6 

100,500 

P 

40 

OH 

J=9/2,F=4-4 

23818.18 

Jan 18, 22 

0.8 

50 ■ 

P 

80. 

OH 

J=9/2,F=4-4 

23818.18 

Jan 22 

0.8 

100,500 

P 

30 

OH 

J-9/2 ,F~5-5 

23826.90 

Jan 23 

0.8 

50 

P 

70 

OH 

- J=9/2,F=5-5 

23S26.90 

Jan 23 

0.8 

100,500 

P 

30 

OH (Comet 

J=9/2,F=5-5 

23826.90 

Jan 23 

0.8 

50 ; i 

P 

40 

Tail) 








CH^OH 

7„-7 '(E) 

25124.88 

Dec 13, 14, 15, 

1.2 

50 

M 

135 


2 1 


Jan 5 




- 

ch 3 oh 

2 5 ,-2 1 (E) 

24934.38 

Dec 12 

1.2 

50 

M 

20 

n 2 0 

1-0 

25123.25 

Dec 13, Jan 5 

1.0 

50 

M 

30 


^Average geocentric distance of comet in A„U. 

2»j — maser preamplifier, P = tunable parametric amplifier 


TABLE XI 


, 

UPPER LIMITS 

FOR MOLECULE 

SEARCHES IN COMET 

KOHOUTEK (1973f ) 

■ -/ 


Molecule 

Transition 

F/R 1 

Peak-to-Peak 

t a (°K) 

Upper Limit 
< T> 

Upper Limit (cm - * 2 ) 
<N> 

H 2° 

6 16 " 5 23 

50 P 

0.80 

.0060 

5.0 x 

10 16 

kh 3 

. (1,1) 

50 M 

0.25 

.0024 

5.4 x 

10 13 

nh 3 

\ 

(2,2) ' 

100 P 

4.07 

.0400 

5.0 x 

10 14 

nh 3 

(3,3) 

50 M 

0,28 

.0027 

1.5 x 

10 13 

OH 

J=9/2,F=4-4 

50 P 

1.35 

,0133 

6.2 x 

10 16 

OH 

J=9/2,F=5-5 

50 P 

0.92 

.0104 

4.2 x 

10 16 

OH (Tail) 

J=9/2,F=5-5 

50 P 

1.75 

.0170 

8.0 x 

• 

10 16 

CH o 0H . 

7 2 - 7 X (E) 

50 M 

0.32 

.0032 

9.0 x 

10 14 

ch 3 oh 

' 2 2 " 2 i( E > 

50 M 

• 0.91 

.0092 

8.4 x 

10 15 

n 2 ° 

1-0 

50 M 

1.89 

.0194 

1.5 x 

10 17 


1 F/H denotes filter resolution in kHz/receiver , where M * maser preamplifier 

and P = tunable parametric* amplifier 


TABLE III 


UPPER LIMITS FOR MEAN COLUMN DENSITY AND MOLECULE PRODUCTION RATE 

COMET KOHOUTEK (1973f) 







Q 



** 


< N 

> 

r 


10 4 km 

r 

:=2 

10 5 km 

Molecule 

o 


o 



h 9 o 

5 x 

10 16 

1 

X 

10 31 

2 

X 

10 30 

NK 3 

1.5 x 

10 13 

4 

X 

io 27 

6 

X 

IO 26 

CH^OH 

9 x 

r-H 

O 

i-l 

1 

X 

io 29 

2 

X 

10 28 

n 2 0 

1.5 x 

17 

1CT 

2 

X 

10 31 

3 

X 

10 30 


< N > = upper limit for mean column density uniformly spread out over antenna 
‘ beam (molecules~cm ) where the number of molecules refers to the 

total number of molecules in all states 
q = upper limit for molecule production rate (molecules-sec ^-steradian 
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ABSTRACT 

An upper limit on the rate of evolution of CO from Comet 
Kohoutek has been obtained from a search for resonant scattering 
of sunlight near The observations were made approximately 

two months after perihelion. The rate of evolution of CO at that 
time was apparently less than that of CH^CN observed before 
perihelion . 



Under favorable conditions, and with a very sensitive 
spectrometer, Infrared resonant scattering of sunlight should 
allow a determination of the rate of evolution of CO from a 
bright comet (Blermann, 1973a). Such a determination was 
attempted for Comet Kohoutek 1973f. Optimum observing condi- 
tions were not obtained and CO not detected, ■ but the limited 
measurements which could be made provide an upper limit on. the 
CO evolution rate approximately two months after perihelion. 

Observations were made on February 23 > 197^ at the Mauna 
Kea Observatory’s 88" telescope using a tandem scanning Fabry- 
Perot interferometer (Holtz, 1971)* The field of view was 9.5 ,T 
by 14” and was approximately centered on the comet’s nucleus. 

A spectral region near 4.7t including the 1-0 P2 and P3 funda- 

12 16 

mental rotation-vibration transitions of CO was scanned with 

a resolution of 0.1 cm” 1 . The comet’s radial velocity of *17 

—1 

km s with respect to the earth, giving a doppler shift of 
0.33 cm" 1 , was sufficient to move the expected positions of the 
cometary lines well outside the half-power points of the terres- 
trial CO line absorption. 

The sensitivity of the spectrometer system was determined 
by measurement "'of the signal from W Hya in the same wavelength 
region. If o represents the r.m.s.- noise fluctuation, any 
signal larger than 3a should have been easily observed. With 
this criterion, the lack of an observed signal sets an upper 



5 


limit of 1.0 x ICT 18 w cm" 2 for the flux In either CO line at 
the top of the earth's atmosphere. Since the rate of resonant 
scattering of 4.7p radiation by a CO molecule at 1.5 A.U. from 
the sun is approximately 6 x 10" 5 s' 1 , the above limit for the 
flux corresponds to an upper limit of 2.5 x 10 33 for the total 
number of CO molecules within the field of view in a single 
rotation-vibration state. For an assumed temperature of 250° K, 
the J= 3 level contains about 1% of all CO molecules present . 
Thus, the upper limit on the total amount of CO in the field 
of view is 3.5 x io 3 i 

Assuming a model in which CO is produced in a region small 
compared to the observed field of view, the density of CO is 
given by n(CO) = -2_ , where ±s the total productlon 

rate, r the distance from the source, and v is the mean speed 

‘ F0r a molecu lar lifetime t, the CO cloud will have a 
radius r Q = fi. The total number of molecules within a circu- 
lar field of radius p centered on the comet is approximately 

of ^TrQ v . 

v P if P<r Q> which is almost certainly the case. At a 
distance of 1.5 A.U., the rectangular field of view has the 
area of a circular field of radius p = 7-3 x 10 8 cm. Hence 
the observed limit on the flux corresponds to Q/v (CO) < 2x10^ 

cm ster 1 . Biermann (1973a) has estimated that Q/v for CO 
in bright comets is 
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Ulich and Conklin (197*0 have reported a detection of 
CH 3 CN in Comet Kohoutek by its microwave emission and concluded 
that the total number of methyl cyanide molecules within their 

■3 /) 

antenna beam was ~9 x 10 . The observations were made in 

early December, before perihelion. They estimate the CH^CN 
cloud radius to be -lO^cm, which is considerably smaller than 
the size of the microwave antenna beam. In that case, 

Q/v (CH CN) = £ = 7.2 x 10 22) cm " 1 ster -1 

^ P 1 j 

K cloud 

If the CH^CN cloud is larger than the estimated lO^crn, a lower 

limit on Q/V is given by -r— . The effective radius of the 

p beam 

beam used by Ulich and Conklin is -4 x 10 y cm, so that 
Q/v (CH^CN) >1,8 x 10 24 cm ” 1 ster ” 1 
This- is essentially equal to the above upper limit for CO, 

Q/v (CO) < 2 x 10 24 cm 1 ster 1 
It is surprising- that the rate of, evolution from the comet 
reported for CH^CN is equal to or greater than the rate of 
evolution of CO. However, it must be noted that CH^CN was de- 
tected before perihelion when the comet was approximately 0.8 
A.U. from the sun. Subsequent search after perihelion failed 
to detect this molecule. Observations which set the present 
upper limit on CO were made only after perihelion, when the 
comet was ~1.5 A.U. from the sun. Since one would expect the 



rate of . evolution of CO to be much, greater than that of CH^CN, 
it is unlikely that this difference in heliocentric distance 
can account for the apparent abundances. It is more likely 
that the comet lost most of its volatile gases such as CO in 
the course of its trajectory before and during perihelion, 

(cf. Biermann, 1973b) or that the more volatile' gases were 
already depleted long before its recent approach to the sun. 
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SYMBOLS 

\i lower case Greek mu 

a lower case Greek sigma 

it lower case Greek pi 

( )" arc seconds (as 9 *5” by 14") or inches (as 88”) 
< less than 

> greater than 

- approximately equal to 

v lower case Arabic v with bar 



